On the effective Dirac dynamics of ultracold atoms in bichromatic optical lattices 
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We study the dynamics of ultracold atoms in tailored bichromatic optical lattices. By tuning the 
lattice parameters, one can readily engineer the band structure and realize a Dirac point, i.e. a true 
crossing of two Bloch bands. The dynamics in the vicinity of such a crossing is described by the 
one-dimensional Dirac equation, which is rigorously shown beyond the tight-binding approximation. 
Within this framework we analyze the effects of an external potential and demonstrate numerically 
that it is possible to demonstrate Klein tunneling with current experimental setups. 

PACS numbers: 03.75.Mn,03.67.Ac,03.65.Pm 



I. INTRODUCTION 

Quantum simulators aim at the simulation of complex 
quantum systems in well controllable laboratory exper- 
iments Q. Such a simulation is especially useful when 
the original quantum system is experimentally not ac- 
cessible and numerical simulations are impossible due to 
the exponential size of the Hilbert space. Furthermore, 
quantum simulators offer the possibility to tune the ex- 
perimental parameters to explore novel physical phenom- 
ena. Important examples include the simulation of solid 
state systems with ultracold atoms sonic black holes 
in Bose-Einstein condensates |3[ and the Dirac dynamics 
with trapped ions 0-0 • 

Ultracold atoms in optical lattices are especially suited 
for such a task, since their dynamics can be controlled 
with an astonishing precision and their dynamics can be 
measured in situ. Bichromatic lattices are especially ap- 
pealing since these systems allow to tune the energy dis- 
persion of the Bloch bands. In particular one can choose 
the parameters such that a Dirac point, i.e. a true cross- 
ing of the first and second excited band, is realized de- 
pending on the relative phase between the two funda- 
mental lattices. Unlike other systems (8l4Toj. bichromatic 
optical lattices thus allow to simulate relativistic quan- 
tum effects using only a single species of neutral atoms 
and no external driving fields. This approach therefore 
paves the way for the simulation of interacting relativistic 
quantum field theories [11] . 

In this paper we investigate the dynamics around a 
Dirac point in detail and derive the one-dimensional 
Dirac equation as an effective equation of motion for the 
coarse-grained atomic wave functions. In contrast to pre- 
vious approaches flll - [l4j , we do not make use of a tight- 
binding approximation, such that the Dirac equation is 
found without imposing a continuum limit. We discuss 
the effects of an external potential in detail, showing that 
the Dirac description remains valid if the potential varies 
slowly enough. Within this framework we finally show 
that it is possible to simulate Klein tunneling through a 
potential barrier with current experimental methods. 



II. BLOCH AND WANNIER STATES IN 
BICHROMATIC OPTICAL LATTICES 

We consider the dynamics of ultracold atoms in a 
bichromatic optical lattice described by the Hamiltonian 



-h 2 d 2 
2M dx 2 



+ — cos(2/cqx) + — cos(4/cqx + qb) (1) 



plus an additional potential V(x), which is assumed to 
vary slowly compared to k$x. A bichromatic optical lat- 
tice with arbitrary relative phase (j) can be implemented 
by a superposition of an ordinary optical lattice with a 
periodicity of A/2 and an additional lattice with a pe- 
riodicity of A/4 based on four-photon processes as de- 
scribed in [l5l-[l7|. For the additional potential V(x) we 
consider (i) an optical dipole trap and (ii) a static field 
which can be realized by either gravity or accelerating 
the complete lattice. In the following we will use scaled 
units which are obtained by setting x' — k§x, t' = Ejit/h 
and dividing the Schrodinger equation by the recoil en- 
ergy Er = h 2 k 2 ,/2M. In these units we have H = 1 and 
M = 1/2 and all energies are given in units of Er. 

The eigenstates of Hq, the Bloch waves can be engi- 
neered to a large extend by choosing the lattice parame- 
ters Vi 5 2 and qb which makes bichromatic lattices a con- 
venient tool for quantum simulations. To be precise, the 
Bloch states are defined as the simultaneous eigenstates 
of Hq and the translation Td over the lattice period d: 



H U a ^(x) = E a (K)Ua iKl {x), 

T d u a , K (x) = e idK u ayK (x). 



(2) 



Here and in the following, k G [— ir/d, ir/d] is the quasi 
momentum and a = 0, 1, 2, . . . labels the different Bloch 
bands. Figure [1] shows the bandstructure E a {n) of a 
bichromatic lattice with V\ = 5 and V2 = 1.56, compar- 
ing two different choices of the relative phase <fi. For these 
values of the lattice depth and a relative phase (j) = tt, 
one observes a true crossing of the eigenenergies of the 
first and second exited band at k = 0. i.e. a so-called 
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FIG. 1: The lowest Bloch bands for the case Vi = 5 and 
V2 = 1.56. For = 7r one observes a true crossing between 
the first and the second excited band, a Dirac point, which 
can be used to simulate the Dirac equation. 



'Dirac point'. The physical reason for the vanishing of 
the band gap is that the contributions of the second or- 
der Bragg scattering at the optical lattice with periodic- 
ity A/2 and the first order Bragg scattering at the lattice 
with periodicity A/4 show a complete destructive inter- 
ference for the given parameters. At the Dirac point the 
dispersion relation is linear in ft, just as for a relativis- 
ts massless Dirac particle, such that a similar dynamics 
can be expected. We will make this analogy more precise 
in the following. In general, the band gap between the 
first and second excited band is approximately given by 
AE « |(Fi/4) 2 + V 2 exp(z0)| 0. 

However, also for a small but finite band gap we obtain 
a pseudo-relativistic dynamics in the center of the Bril- 
louin zone. In any case, we can approximate the energy 
dispersion of the first and second excited band around 
ft « as 



£i,2(V) = e d ± \/m 2 c 4 



(3) 



This relation defines an effective mass m which is given 
by the curvature of the two bands and an effective speed 
of light c which is related to the band gap by 

AE = 2mc 2 . (4) 

The applicability of this approximation is illustrated in 
Fig. [2j where it is compared to the numerically exact 
data. Furthermore, the figure shows the effective param- 
eters m and c as a function of the lattice phase </>. 

Now it is very convenient to introduce a new basis in 
which the two Bloch states in the first and second excited 
band are rotated: 



U2,k = - smQui^ + cos# u 2 ,k 
with the mixing angle 



tan#(ft) 



rnc 



; + Vm 2 c 4 + c 2 ft 2 



(5) 
(6) 
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FIG. 2: (Color online) (a,b) A fit of the relativistic dispersion 
relation Q to the Bloch bands a = 1 , 2 in the center of the 
Brillouin zone for cj) = and cj) = 7r, respectively. (c,d) The 
resulting fit values for the effective parameters m and c as a 
function of the phase <j). The remaining parameters are the 
same as in Fig. [T] 



In this basis, the free Hamiltonian i^o is no longer diag- 
onal, but has the convenient form 



r r , x ( E D + ck rnc 2 
H ^)=\ mc 2 E d -ck 



The eigenenergies (|3]) for given quasimomentum k are 
simply the eigenvalues of the matrix Hq(k). 

For an effective description of the quantum dynamics in 
a perturbed crystal we will furthermore need the Wannier 
basis which is defined as follows. The Bloch waves can 
be chosen to be periodic in the quasi momentum ft, such 
that they can be expanded into a Fourier series, 



u aiK (x) = -^Ve iK(i \ )n (x), 



(8) 



which defines the Wannier states w^ n {x). Inverting the 
Fourier series yields 



w ay7l (x) 



— I 

V2^d J 



e~ indK u a ^(x) dn. (9) 



The Wannier states are exponentially localized at the lat- 
tice site n [19] , and they are related by a simple shift in 
real space 



(10) 



where x n = nd is the position of the nth lattice well. 
These properties are quite useful for several approxima- 
tion schemes (cf. [20]). Examples for the case of a bichro- 
matic lattice are shown in Fig. [3] for two values of the 
relative phase <p. 



3 




2 

x/d 



FIG. 3: Squared modulus of the Wannier states \w a , n (x)\ 2 
in the ground (a = 0, upper panels) and the first excited 
band (a = 1, lower panels) at the lattice site n — in a 
semi-logarithmic scale. We assume a bichromatic optical lat- 
tice with Vi = 5, Vb = 1-56 and = (left) and cj) — 7r, 
respectively. 



III. EFFECTIVE EVOLUTION EQUATIONS 

For the derivation of the effective evolution equations 
we start from the Wannier representation of an arbitrary 
wave function, 



^(x) = ^ ^<x,nW a ,o(x - X n ) 



(11) 



The expansion coefficients t/j a ,n form an infinite, but 
countable set of complex number. For the following ap- 
proximations, however, we define a continuous function 
ipatx) for every Bloch band a such that 



(12) 



These functions can be viewed as a coarse grained ver- 
sion of the quantum state ^(x) projected on the band 
a. Now one can derive effective evolution equations for 
these coarse grained wave functions. 

One can show that every operator, which is diagonal in 
the quasi momentum ft, can be expressed in a very con- 
venient way for the coarse grained wave functions ^ a (x) . 
So consider an operator which satisfies 



(13) 



where u ayK are the Bloch states defined above. Then this 
operator acts on the coarse grained wave function as 

6^(x) = J20eAp)M*)i (14) 

i.e. the quasi momentum n is replaced by the momentum 
operator p = —id x . To carry out this replacement one 
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FIG. 4: (Color online) Matrix elements V a ,n;i3,m of a linear 
potential V(x) — x in a Wannier basis as a function of the 
lattice depth Vi. We have plotted the matrix elements for 
a = 1, P = 2 and \n — m\ < 2, which are the leading cor- 
rections to the diagonal approximation ([19ft . As above, we 
assume V2/V1 = 1.56/5 and = (left) and <f> = tt (right), 
resepectively. 



can expand 0/3^(k) in a Taylor series and then replace 
every term n n by (— id x ) n . This relation was first shown 
by Slater [22| , cf. also |23l-[25j . The proof is summarized 
in the appendix. 

In particular, this holds for the lattice Hamiltonian Hq, 
which is trivially diagonal in the Bloch basis. Going to 
the rotated basis defined in equation (|5j), the eigenvalue 
equation for the lattice Hamiltonian reads 



Ejj + ck rnc 
mc 2 E 



d ~ cn 



(15) 



The constant energy offset Ed introduces a global phase 
shift only, which has no physical significance. By shifting 
the energy scale, we can set it to zero, Ed = 0. 

Now if we express a general quantum state in the first 
and second excited band as 



^( X ) = i^l^a.nWj 



(16) 



where the rotated Wannier function w a ,n( x ) are related 
to the original ones w a ^ n (x) by the same rotation as the 
used in Eqn. (j5j), then Slater's theorem tells us that the 
lattice Hamiltonian acts onto the coarse grained wave 
functions ip a (x) as 



i/j 2 (x) 



r+cp 

Kmc 



rnc 
—cp 



^i(x) 
i/j 2 (x) 



(17) 



where p = — id x is the momentum operator. 

In addition we need to know how a slowly varying po- 
tential V(x) acts onto the functions ip a (x). We recall 
that these functions have been defined as expansion co- 
efficients in the Wannier basis, such that we need the 
matrix elements 



V(3,n;a>,m = J dx Wp {x - X m )V (x)w^ (x - X n ) . (18) 

Using the strong localization of the Wannier states, one 
can set V(x) to a constant over the localization length 
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in a first approximation. Because of the orthogonality of 
the Wannier states one thus finds 

V(3 jn ;a,m ~ V \x n )8 (3 ,a^n,m- (19) 

Let us analyze this approximation in more detail for 
the case of a linear potential V(x) = Fx. Then we have 

Vi3,nia,m = FXn5(3,aSn,rri + F J dx W^ m _ n (x) X W a ,o(x), 

where Eqn. ([TO]) has been used to simplify the results. 
The first term in this expression corresponds to the di- 
agonal approximation ([T9]h while the remaining matrix 
elements can be interpreted as follows. The local terms, 
i.e. the terms with n = m, vanish exactly for a = f3 due 
to the parity of the Wannier functions. The term a = 1 
and /3 = 2 is the most important correction to the diag- 
onal approximation ([T9]) . It couples the two bands and 
can thus be viewed as an additional contribution to the 
effective mass mc 2 in the effective wave equation ([T7]) . 
The nonlocal terms n ^ m are also largest for a = 1 and 
/3 = 2. However, they vanish exponentially with the lat- 
tice depth Vi,2 and the squared lattice period d 2 . These 
matrix elements are plotted as a function of the lattice 
depth in Fig. [4j assuming V2/V1 = 1.56/5 and <j) = 0,7r 
as above. Note that the lattice depth in scaled units is 
proportional to d 2 . One clearly sees how the non-local 
terms vanish exponentially in contrast to the local term 
n = m. 

In the following we will confine ourselves to the first 
order approximation ([T9]) . For the ansatz ([T6]h the po- 
tential thus acts as 

V(x) ^2^a(x n )w ay n(x) ~ ^ V(x n )^ a (x n )w ayn (x) . 

In this approximation, we thus find the effective evolution 
equations for the coarse grained wave function: 

d_ AM _ (V{x)+cp mc 2 \ AM ( . 
dt\^2j~\ mc 2 V(x) -cp)\^ 2 )' [Z{J) 

If we rotate the 'spinor' wave function (^1^2) once again 
by the unitary transformation 

we finally obtain the Dirac equation in 1+1 dimensions 
with an external scalar potential, 

■ d_ (ip a \ _ (V(x)-mc 2 cp \ (ip a \ , 
dt\AJ~\ cp V(x) + mc 2 ) \ih J ■ {ZZ > 

Note that in this rotated frame ?p a and ipb coincide with 
the amplitudes in the first and second excited band for 
ft = - but only there. Otherwise one has to be very 
careful when interpreting the wave functions. 



IV. QUANTUM SIMULATION OF KLEIN 
TUNNELING 

As an example of the effective relativistic dynamics we 
consider the tunneling of a wavepacket out of a dipole 
trap in a tilted bichromatic lattice. In particular, the 
slowly varying potential is given by 

V(x) = -V exp(-2x 2 /^o 2 ) - Fx ( 2 3) 

with F = 0.076, V = 19.77 and W = 157 in scaled 
units, which corresponds to typical experimental param- 
eters [IB, 03] • Initially, the wave packet is localized in the 
second exited band with quasi momentum k = 0.95 with 
a gaussian envelope of width a = 17. 

The resulting dynamics of the atoms is shown in Fig. [5] 
The sketch on the left demonstrates the two different 
regimes realized for <p = (al, bl) and <p = tt (a3, b3). 
In the latter case a Dirac point emerges in the band struc- 
ture, such that the atoms behave like massless relativistic 
particles. Unlike a massive Schrodinger particle such a 
Dirac particle can escape from the trap via Klein tunnel- 
ing. 

This expectation is confirmed by the numerical simu- 
lation of the atomic dynamics. The upper panels (al- 
a3) of Fig. [5] show the evolution of the modulus of the 
atomic wave function calculated with the origi- 

nal Schrodinger equation in a bichromatic optical lattice. 
The effective Dirac dynamics according to Eqn. ([20]) is 
shown in the lower panels (bl-b3). The effective val- 
ues of the mass m and the speed of light c are given by 
mc 2 = 0.78 (bl), mc 2 = 0.24 (b2), and mc 2 = (b3), 
respectively. For the sake of a better visibility we have 
again plotted the modulus \ipi(x, t) \ + \?p2(%,t)\. One ob- 
serves that the essential features of the atomic dynamics 
are very well reproduced by the Dirac approximation. 
In particular, one observes the transition from a 'heavy' 
Schrodinger-like particle for <p = (al, bl) to a 'rela- 
tivistic' particle with vanishing effective mass for </> = tt 
(a3, b3), which escapes from the dipole trap by Klein 
tunneling. For <j> = 0.8 7r (a2, b2), partial tunneling is 
observed and the potential barrier acts as a matter wave 
beam splitter. 

The most obvious difference of the Dirac approxima- 
tion ([20]) to the underlying lattice dynamics is that the 
group velocity of the Dirac wave packet is limited to the 
effective speed of light c. After tunneling out of the dipole 
trap, the atoms are accelerated by the linear potential 
Fx, which is not observed in an effective relativistic de- 
scription. This difference is due to the fact that the wave 
packet is not restricted to the center of the Brillouin zone 
for the given parameters, as it was assumed in the deriva- 
tion of the effective Dirac equation. Instead, we have 
^initial = 0-9 and |ftfinai| < 1 in the example shown in 
Fig. [51 
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FIG. 5: (Color online) Quantum simulation of Klein- Tunneling out of a dipole trap. The upper panels (al-a3) show the full 
time evolution of the Schrodinger equation in a bichromatic optical lattice for the case Vi = 5 and V2 = 1.56 and a relative 
phase of cj) = (al), 4> = 0.8 ty (a2), and (f> = ty (a3). The lower panels (bl-b3) show the dynamics of the effective Dirac equation 



([20|) with effective mass and speed of light given by 



0.78 (bl), mc 2 =0.24 (b2), and mc 2 = (b3), respectively. 



CONCLUSION AND OUTLOOK 
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We have analyzed the quantum dynamics of ultracold 
atoms in a bichromatic optical lattice. It was shown that 
the lattice parameters can be tuned such that a Dirac 
point emerges in the band structure, i.e. a true crossing 
of the Bloch bands with linear dispersion relation. In the 
vicinity of such a crossing the atoms effectively behave 
like massless relativist ic particles, allowing for a tabletop 
simulation of relativist ic quantum physics. 

We have rigorously shown that the one-dimensional 
Dirac equation is found as an effective evolution equation 
for the coarse grained atomic wave function projected 
onto the two crossing Bloch bands. Unlike previous ap- 
proaches, our derivation does not rely on a tight-binding 
approximation and also shows that how to include an 
additional, slowly varying potential. Therefore it is pos- 
sible to simulate Klein tunneling - the tunneling of an 
ultrarelativistic particle though a potential barrier with- 
out damping - with current experimental setups. 

Ultracold atoms in bichromatic optical lattices have 
some important advantages than other systems proposed 
before. The initial state of the atoms and all experimen- 
tal parameters can be controlled with astonishing preci- 
sion. In a common experiment, several thousands of ul- 
tracold atoms are prepared in the optical lattice. Thus it 
is possible to simulate the dynamics of interacting Dirac 
fermions in contrast to experiments with single trapped 
ions. While these experiments of course require the use 
of fermionic ultracold atoms, we note that Klein tunnel- 
ing is a single particle effect such that it can be observed 
equally well for bosonic atoms. 
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Appendix: Slater's derivation of the effective 
equation of motion 

The wave function ^(x,t) is expanded into Wannier 
states 

= ^ j>ot (Xn)Wa,0 (x - X n ). (A.l) 



The wave function is thus represented by a discrete set of 
numbers ^ a ,n- Due to the exponential localization of the 
Wannier states, these coefficients can be interpreted as 
the amplitude in the nth lattice site and in band a. For 
the effective evolution equations treat these coefficients 
as a continuous function in x, i.e. we choose a smooth 
function ip a (x) such that 



Safari) — ^a,n- 



(A.2) 



This function can be seen as a coarse grained version of 
the original wave function ty(x) projected onto the ath 
Bloch band. 

We consider an operator, which is diagonal in the quasi 



6 



momentum 

6 U a , K = ^ °/3,aW^ )/t , (A.3) 

/3 

as for instance the unperturbed Hamiltonian Ho. The 
functions Op, a (K,) are periodic in ft, such that we can 
expand them into a Fourier series: 

0^) = Y. 0i pl e ~ iSdK - (A- 4 ) 

S 

Applying this operator to a wave function of the form 
(jA.l|) and inserting the definition ((9]) for the Wannier 
function yields yields 

(x) 




Using again equation (|9]) this can be rewritten as 

6* (a?) = ^ ^a(Xn)0$ a Wp(x - X n + 8 ) 
ac,/3,n,s 

= ^2 °$a^oc{Xrn- Xs)™p{x - X m ), 
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